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Abstract 

Structural information in the local electric field produced by helical charge distribu- 
tions, such as dissolved DNA, is revealed in a straightforward manner employing cylin- 
drical coordinates. Comparison of structure factors derived in terms of cylindrical and 
helical coordinates is made. A simple coordinate transformation serves to relate the Green 
function in cylindrical and helical coordinates. We also compare the electric field on the 
central axis of a single helix as calculated in both systems. 
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A few years ago we developed an exact analytical solution for a model of the local electric 
potential and field arising from the double helix of phosphate groups of a single B-DNA molecule 
immersed in an aqueous solvent M . We subsequently extended our calculation to treat the full 
(sugar-phosphate plus base-pairs) discrete charge distribution of homopolymer B-DNA in a 
solvent modelled by concentric dielectric cylinders [Q. In both cases we found a characteristic 
length scale of ~ 5A associated with the radial persistence length of either the helical imprint 
or individual base-pair identity. These detailed calculations are based on a theoretical model 
of B-DNA in solution together with Green function techniques to account for the contribution 
of each individual (partial) charge to the full net potential. As we discussed in detail in g], 
although the helical configuration of point charges does not possess cylindrical symmetry, the 
individual charges making up the backbone as well as the base-pairs can be assigned to a set of 
concentric cylindrical surfaces and thus cylindrical coordinates provide the most natural system 
in which to separate Laplace's and Poisson's equations. The full power of the Green function 
technique can be brought to bear on the problem, which is made nontrivial by the presence 
of a non-uniform dielectric. Moreover, in cylindrical coordinates one may easily identify and 
separate out the featureless zero-mode contribution to the potential. The zero mode, which 
goes as ~ ln(p), corresponds to the cylindrically symmetric potential from a straight line of 
charge and is what one sees far from the DNA surface. The higher-mode terms in the potential 
therefore encode and reveal the specific helical conformation of the molecule. The structural 
information contained in these higher modes is complete, transparent and calculable [jl], |2[ . 

Though the above points have been adequately discussed in |], §, we wish to emphasize 
them in light of a recent claim that structural information of (double) helix charge distributions 
is better represented employing helical rather than cylindrical coordinates ||. There, it was 
argued that structure factors derived in a certain helical coordinate system reveal structural 
information in a more "transparent" fashion. The purpose of the present brief report is to 
demonstrate that similar structural information may be represented just as easily in cylindrical 
coordinates. 

Structure Factors 

We turn to the derivation of structure factors in cylindrical coordinates. In the interests of 
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clarity and brevity, we shall restrict our attention to the case of a single helix immersed in a 
single dielectric background medium. In our two previous papers we calculated explicitly both 
the electrostatic potential and electric field components due to many helices (to account for 
the phosphate backbones and base pairs) immersed in a medium modeled by a triply piecewise 
constant dielectric function. However, since the single helix in a uniform background is the only 
case treated by the methods in [[| , we will stick to this simple example to facilitate comparison. 

Before tackling the full helix, we begin by considering the electrostatic potential for a single 
point charge q embedded in a uniform dielectric medium e, and located at the point (a, </>', z'). 
^From 0] or the expression (B26) in the Appendix of [I], we have 



where I m} K m are the modified Bessel functions of integer order m and the prime on the sum 
indicates that the zero mode is to be divided by two. ^From this we can build up the electro- 
static potential corresponding to a finite or an infinite distribution of (identical) point charges 
distributed along a helix of constant radius a and pitch P. We regard a single helix of evenly 
spaced point charges as a two-dimensional regular lattice wrapped around the cylinder p = a. 
Instead of summing the charges along the helix, we decompose this equivalent lattice into a 
finite collection of one-dimensional vertical line charges that run parallel to the z-axis and then 
sum over this "bundle" of one-dimensional chains. Denote the number of such chains by the 
integer N Q and the number of charges along any such chain (this can be finite or infinite) by 
2N + 1. The total number of charges living on the helix is therefore equal to N x (2N + 1). 
We now make this conceptual decomposition quantitative: the position in cylindrical coordi- 
nates of the n-th charge on the s-th vertical line with fixed angular coordinate 4> s is as follows: 
(a,<p s ,Zn iS ), where 



z n)S = nP + sAz, 

for < s < N — 1 and — N < n < N, where Az is the vertical rise per residue, and N Q = P/Az 
is the number of equally spaced residues per pitch length of the helix (see Fig. 1 of ref. [1]). 



$(p,0,^) = ^E / 



•oo 



dk J m (/cp < )i^ m (/cp > ) cos(k[z — z'))) cos(m[0 — (/)']) 





(2) 
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The single helix potential is obtained simply by replacing 

z' — > z niS , and d>' — > c 



in (1) and summing over the point charges: 

N No-l 

$Heiix(p,<fr,z) = ^(p,<P,z\a,(f) s ,z njS ). (3) 

n=—N s=0 

This leads to the following expression for the electrostatic potential from a single charged 
helix in a uniform background: 

2q 

"m=0 ,/u ± 

+ S ± (m, k) sin(kz ± m0)}/ m (/cp < )ii' m (/cp > ). 



2q 00 ' r°° 

§Heiix(p,0,z) = — / dkG 1 (k)J2{C ± (m,k)cos(kz±m ( j)) (4) 



The momentum dependent structure factors appearing here are defined as follows: 

= M(2N+im orJvflnite) 

sm(^-J 

00 



and, 



% t 6(k-^f), (foriV^oo), 
i=-oo 
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C^im^k) = Gf{m,k)cos[{N -l)(k±—p-)Az/2), (6) 

2 7T/77 

S ± (m,k) = G$(m,k)em[(N -l)(h±-j r )Az/2]. 

The remaining structure factor appearing in (4) is defined by 

sin[|(LP ± 27rm)] 
G2(m ' fc) = sin[^-(fcP± 27 rm)] - (?) 

The final form in which we have written (4) results from working out the double sum in (3) 
and using one or more elementary series and/or trigonometric identities, as well as the Poisson 
summation formula, as for example, in (A3) of [|TJ. 

By inspection of (4), it is clear that structural information (as defined by the appearance 
of these so-called structure factors) of helical charge distributions is revealed in the cylindrical 
coordinate system. The structure factors depend explicitly on the pitch P and the rise per 
charge, Az, which jointly parametrize the geometry of the helix. 



Note that our G\(k) has the same mathematical form as the function H 2 (M) calculated 
in ||. The reason for this coincidence is clear: both structure factors result from performing 
the charge sum (whether finite or not) over (half) the coordinate lines corresponding to each 
coordinate system. That is, we obtain G\ from summing charges along the coordinate line 
6 = const, (a line parallel to the z-axis) whereas Hi results by summing along the coordinate 
line t = const., that is, along a helix. The linear momentum in the z-direction is k, while that 
in the helical direction t is M. One has in fact the formal correspondence between summing 
over vertical lines and summing over helices: 

k «-> M, (8) 
P <-> s. 

In the helical coordinate system, H2(M) is the full structure factor for a single helix. In 
cylindrical coordinates, G\(k) gives the structure factor for a line charge, and is not the complete 
"answer", as it were. That is why we need the additional factors C (m, k) and S ± (m,k), as 
can be checked, by summing over the N Q vertical chains. In the cylindrical system, the helical 
structure is encoded in terms of the factorized set of structure factors Gi(k)C ± (m, k) and 
Gi(k)S ± (m, k). The necessity of the C ± ,S ± could have been anticipated on general grounds 
if one thinks of Fourier analysis on a cylinder, namely, one expands an arbitrary function of 
and z in terms of the complete set of functions of chiral "up" and "down" moving "waves": 
{sin(A;z ± m0), cosikz ± m(f))}. These coefficient functions (C ± , S ) are simply the projections 
of helical structure along the two helical directions, but expressed in cylinder coordinates. This 
is why these basis functions appear in our $ Helix- We had already made use of this fact during 
an intermediate stage of our published calculation (see, e.g., Eq. (B8) in ffl). 

We may derive another form for the structure factor by performing the sum over charges in 
a distinct way. For example, if we now sum the point charges along the helix, then an equally 
transparent factor results-even using cylindrical coordinates. Only a single sum need now be 
carried out because the helix sum means we are to make the replacements 

0' -> n = nA0, (9) 
z' — > z n = nAz, 
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in (1) for — M < n < M (not be be confused with the momentum variable M appearing in 
(8)) and sum over n to obtain the corresponding single helix potential. The origin charge 
corresponds to n = and has coordinates (a, 0, 0). Note that for a finite helix, the net charge 
must turn out to be the same no matter how we carry out the sum. This implies that 2M + 1 = 
N Q x (2N + 1). By going through identical steps as above, we find that the single helix potential 
is now represented as 

®Heiix(p, 0, z) = — V / dk V{G^(m, k) cos(kz ± m0)}/ m (A;p < )ir m (/cp > ). (10) 
71-6 m=o Jo ± 

Here, the structure factors turn out to be as follows: 

3V ' sin[l(kAz ± mA<p)] K ' 

for finite M (and hence, finite N) or, 

oo 

Gf(m,k) = 2n 5 ( kAz ± mA( P ~ 27l j), (12) 

j=-oo 

for M (or N) infinite. Again, the basic parameters Az, P describing the helix geometry are 
manifest in the cylinder-based structure factors. Moreover, the cylinder and helix-based struc- 
ture factors G3 and H 2 , are mathematically identical in form. The formal replacement is 

Ms -> kAz±mA(f), (13) 

and allows us to obtain one from the other. 

The extension to double helices is straightforward. One merely adds in the potential for 
the dyadically related charges and then sums over the jth backbone charge, q\ making up the 
group: 

®double-hdix(P, <t>, Z ) = ®helix(P, (f> + Oij , Z + Sj) + & helix (p, <f> - <Xj,Z - 6j), (14) 

where 2a j and 28 j are the offset angles and distances, respectively [l], . 

Relation between the helical and cylindrical Green functions 

^From the cylindrical Green function, that is, the potential for a unit (q = 1) point 
charge in vacuum (e = 1) is given by 

-| n oo poo 

- £ / dke- m ^cos[k(z-z')}I m (kp<)K m (kp>). (15) 
_ „ Jo 



* m=-oo J0 
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A similar representation for the left hand side of (15) was given in || which involves the same set 
of basis functions (essentially products of exponentials and modified Bessel functions). Thus, 
a coordinate transformation should suffice to relate these two Green functions. Since we are 
talking about the same function, and provided both coordinate representations are correctly 
worked out, then it must be the case that 

1 , 1 



\ cylindrical i ,i \ helical i 



(16) 



to be understood as an equality between one and the same function expressed in two different 
coordinate systems. Extending the range in the variable k in the above Green function (15) by 
replacing k —>■ \k\ in the arguments of the Bessel functions, use 

/ dk cos[k(z - z')} = - dke- ik{z - z \ (17) 

to rewrite the above Green function (15) (see also Q). Now make (or define) the coordinate 
transformation (0, z) — * (t, s): 

z = s sin (3 — t cos (3, (18) 

s t 
<f) = - cos j3 H — sin p, 
a a 

which is just the transformation from cylindrical to helical coordinates defined in ||. Here, f3 
is the pitch angle of the helix, and is related to the helix pitch by P = 27ratan/3. This puts 

provided we identify 

in 

M = — cos p + k sin /3, (20) 
a 
vn 

M' = — sin P - k cos p. 
a 

This identification can itself be inverted to yield 

m = a{M' sin/3 + Mcos/5), (21) 
k = -M'cosP + Msinp. 

In (21), m is an integer while k represents the linear momentum along the ^-direction. These 
are written in terms of M' and M, which represent the components of linear momentum along 



the helical t and s-directions, respectively. Using (17), (19) and (21), we can now rewrite the 
helical representation for r^z^/i calculated in || and transform it back in terms of cylindrical 
variables. That expression was given as 

, 1 , = - r dMdM'e-^-^e-M'^hikp^Kxikp^, (22) 

|X — X'| 71 J-oo 

where A = a(M' sin (3+ M cos (3) and k is as above in (21) . The first step is to insert the identity 

/oo 
dm 5 (m - aM' sin (3 - aM cos/?) = 1, (23) 
-oo 

and perform the integration over M'. This gives a constant factor (a sin/?) -1 ; this step is 
followed by replacing the integration over m by a discrete sum, / — — ► Yl m i i n accord with the 
fact that m/a is the component of linear momentum in the 0-direction of the cylinder. Next, 
the identity 

f°° j, r/i aM -mcos(3. 

/ dkS{k . t, ) = 1, 24 

j-oo a sm p 

allows one to carry out the integration over M. This produces a factor of sin f3 which cancels 
the (sin/3) -1 from the previous step. Carrying out these steps together with the relations in 
(19) and (20), demonstrates the equivalence (16) between the two representations calculated 
for the Green function. 
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Electric field on the central axis 

Another point of direct comparison with the results of || is had by considering the calcula- 
tion of the electric field on the central axis of the helix. It is not necessary to employ complicated 
Green functions in order to calculate the components of the electric field on the central axis 
in a uniform medium. The potential for a single point charge located at x' = (a, <//, z') in a 
uniform dielectric medium is just 

= iwh-y (25) 

q 1 



e (p 2 + a 2 - lap cos(</> - 0') + (z - z 1 ) 2 ) 1 / 2 ' 
A textbook calculation yields the electric field on the central axis: 



g„ acos(^ 27rsAz * 



E ^' Z > = -^ (,H (: -,F-^y (26) 



(a 2 + (z-nP-sA 2 ) 2 ) 3 / 2; 

(c\ a \ ?r (z-nP - sAz) 
E z (0,(j),z) = -2^ 



e^(a 2 + (z-nP-sAz) 2 ) 3 / 2 ' 

Since these formulas involve the features of the structure of the DNA (clearly through the 
two parameters and P), their dependence on the lengths and nature of the DNA strands 
are also obtained here explicitly. 

The much lengthier calculations of the electric field components on the (double)-helix central 
axis in || lead to the same results only after that calculation is corrected for a spurious factor 
of sin/3 0. Indeed, the step invoking the replacement of integration over M' in (22) by a sum 
over A(= m) in |§ overlooked this factor of (sin/?) -1 which results, as we see, from correct use 
of the identity (23). Once this is done, the two sets of electric field components as calculated 
in terms of cylindrical and helical coordinates do agree. 

Some concluding remarks are in order. First of all, and most importantly as demonstrated 
above, not only are cylindrical coordinates adequate for calculating potentials and fields due to 
helical charge distributions, they also lead to simple structure factors admitting a straightfor- 
ward physical interpretation. 
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In regards to the helical coordinate system advocated in ||, it is a fact that Laplace's 
equation does not separate in helical coordinates. Indeed, it has been known for some time 
that there are exactly eleven coordinate systems in which the three-dimensional Laplacian is 
separable ||; the helical system is not one of them. Nevertheless, the helical Laplacian in |3j was 
separated assuming an a-priori, but not general, form for the t and s-dependent factors of the 
product ansatz. For this reason, it is important to establish the validity of the Green function 
so constructed and the coordinate transformation between cylindrical and helical systems leads, 
as shown above, to independent confirmation of the correctness of the helical Green function. 
At this stage, in so far as one restricts attention to charged helices in uniform dielectric media, 
which coordinate system one prefers to employ is simply a matter of taste. However, for the 
more realistic case of a non- uniform dielectric or solvent, the cylindrical coordinate system offers 
an advantage |TJ, 0. In this situation the question of the matching of solutions at the dielectric 
interfaces arises. Charges within a dielectric induce a surface charge at the interface between 
that dielectric and a region with a different dielectric constant. Since helical coordinates are 
non-orthogonal and curvilinear we expect that matching solutions of Laplace's or Poisson's 
equations across the boundary will lead to a nontrivial albeit solvable problem in differential 
geometry. We in fact had already considered, but rejected, using helical coordinates in 1993 for 
some of the above reasons. Finally, even though helical coordinates are not separable, we can 
always take the potential and electric field computed in cylindrical coordinates and then simply 
transform them to helical coordinates to find their correct forms in these latter coordinates. 
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